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ABSTRACT 

We propose a strategy to search for bulk motions in the intracluster medium (ICM) of merging clusters based 
on Chandra CCD data. Our goal is to derive robust measurements of the average redshift of projected ICM 
regions obtained from the centroid of the Ka line emission. We thoroughly explore the effect of the unknown 
temperature structure along the line of sight to accurately evaluate the systematic uncertainties on the ICM 
redshift. We apply our method to the “Bullet cluster” (IE 0657-56). We directly identify 23 independent 
regions on the basis of the surface brightness contours, and measure the redshift of the ICM averaged along 
the line of sight in each. We find that the redshift distribution across these regions is marginally inconsistent 
with the null hypothesis of a constant redshift or no bulk motion in the ICM, at a confidence level of about 
2 cr. We tentatively identify the regions most likely affected by bulk motions and find a maximum velocity 
gradient of about (46 +13) km s ' kpc ' along the line of sight on a scale of ~ 260 kpc along the path of the 
“bullet”. We interpret this as the possible signature of a significant mass of ICM pushed away along a direction 
perpendicular to the merging. This preliminary result is promising for a systematic search for bulk motions 
in bright, moderate-redshift clusters based on spatially resolved spectral analysis of Chandra CCD data. This 
preliminary result is promising for a systematic search for bulk motions in bright, moderate-redshift clusters 
based on spatially resolved spectral analysis of Chandra CCD data. 

Subject headings: galaxies: clusters: intracluster medium — galaxy clusters: individual (IE 0657-56) — X- 
rays: galaxies: clusters 


1. INTRODUCTION 

X-ray spectral diagnostics of galaxy clusters is a very pow¬ 
erful tool for characterizing the properties of the intraclus¬ 
ter medium (ICM). In addition to density, temperature, and 
metal abundance of the diffuse baryons, another piece of in¬ 
formation is the velocity derived from the Doppler shift of the 
spectral lines emitted by heavily ionized elements. ICM ve¬ 
locity is a critical diagnostic to investigate the dynamics of 
the cluster and its evolution. In particular, bulk motions and 
turbulence are expected as a consequence of major merger 
events. The majority of the merger events occurring during 
the lifetime of a galaxy cluster can produce gas motion at the 
level of 100 km s ' (Nagai et al. 2013). However, major, off- 
center merger events are rare but may provide significant an¬ 
gular momentum to the ICM associated with bulk velocities 
of several thousand km s '. 

The ubiquitous X-shell line complex of the H-like and He- 
like iron in the 6.7-6.9 keV rest frame, originally detected by 
Mitchell et al. (1976), is the most prominent feature in the 
X-ray spectra of massive clusters. Many other emission lines, 
including the L-shell from iron and a element transitions, may 
be found particularly in the soft energy range (0.5-2.0 keV). 
However, considering that clusters are now detected in the X- 
ray up to redshift z ~ 1.75 (see Brodwin et al. 2015), the Ka 
iron complex is the only one which is potentially detectable 
in the ICM at any redshift with a total number of net counts 
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as low as ~ 1000 in the entire X-ray band. Indeed, the 6.7- 
6.9 keV emission-line complex has been identified in the large 
majority of the clusters observed in X-ray at z < 1.5 (Rosati 
et al. 2004; Stanford et al. 2005; Rosati et al. 2009; Tozzi et al. 
2013), while this appears to be challenging at z > 1.5 (see, 
e.g., Tozzi et al. 2015). 

Needless to say, the identification of iron line emission is 
crucial for measuring the redshift of newly discovered clus¬ 
ters in X-ray surveys. Yu et al. (2011) searched for the X„- 
shell iron emission lines with blind X-ray spectroscopy in 46 
Chandra clusters in order to investigate the reliability of the 
measurement of X-ray redshift in medium- and high-z clus¬ 
ters. The goal was to find the data quality requirements nec¬ 
essary to achieve robust redshift measurements in future X- 
ray cluster surveys. Recently, Tozzi et al. (2014) applied this 
technique to measure the redshift of a few newly discovered 
clusters in the Swift X-ray Cluster Survey (Tundo et al. 2012; 
Liu et al. 2015). 

As already mentioned, the Doppler shift of the X^-shell 
emission lines of iron also offers a direct means of investigat¬ 
ing the dynamics of the ICM along the line of sight through 
spatially resolved X-ray spectroscopy. Despite this interest¬ 
ing possibility, the iron line complex has rarely been used to 
search for bulk motions in the ICM. The main reason behind 
this is that the spectral resolution of CCD spectra at the rele¬ 
vant energies is too low compared to the expected effects. In 
the case of Chandra, the EWHM is ~ 285 and ~ 150 eV in 
ACIS-I and ACIS-S, respectively, at 5.9 keV Considering 
only low-redshift clusters with z < 0.3, the resolving power 
in the observing frame around the iron line complex is in the 
range 18-24 and 35-46 in ACIS-I and ACIS-S, respectively. 
This allows one to measure velocity differences of the order 
of at least a few thousand km s ', which may occur only in 

^ http://cxc.harvai'd.edu/proposer/POG/htmI/chap6.html 


2 


major mergers. On the other hand, grating data with much 
higher spectral resolution power, of the order of ~ 1000, can 
only probe the soft band, and will only have angular reso¬ 
lution in the cross-dispersion direction, making it difficult to 
identify regions with different motions. 

Despite these difficulties, CCD data has the advantage of 
providing spatially resolved spectral analysis, which is cru¬ 
cial to unambiguously separate ICM regions with different 
velocities. A few cases of successful measurements of bulk 
motions in the ICM with CCD data have been reported in the 
literature. A first attempt concerned the Centaurus cluster, 
where bulk motions in the ICM were found in the spatially re¬ 
solved spectral analysis of ASCA (Dupke & Bregman 2001a) 
and, eventually, of Chandra data (Dupke & Bregman 2006), 
for a velocity difference of (2.4 + 0.1) x 10^ km (but see 
also Ota et al. 2007). Signs of bulk motions were also found 
by ASCA in the Perseus cluster (Dupke & Bregman 2001b). 
A systematic search with ASCA found reliable signatures of 
bulk motions in the ICM in 2 out of 12 low-redshift (z <0.13) 
clusters (Dupke & Bregman 2005). A significant difference in 
velocity of (5.9 + 1.6) x 10^ km i"' between two regions of 
Abell 576 has been found by combining Chandra and XMM- 
Newton data (Dupke et al. 2007). Other clusters have been 
simply reported to show a significant discrepancy between the 
average redshift of the ICM and that of the BCG, as in the case 
of Abell 85 observed with XMM-Newton (Durret et al. 2005). 
This can be interpreted as a signature of different dynamics 
between the collisionless stellar component and the ICM, and 
it does not represent a measurement of bulk motions within 
the ICM. 

A few single clusters were studied in detail with Suzaku 
thanks to the spectral resolution of the X-ray Imaging Spec¬ 
trometer, which is comparable to Chandra ACIS-I. However, 
Suzaku has a better sensitivity to the Ka iron emission-line 
complex than Chandra and XMM-Newton thanks to its lower 
background. Only upper limits of the order of 2000 km i ' 
are derived for AWM7, Abell 2319 and Coma (Sato et al. 
2008; Sugawara et al. 2009; Sato et al. 2011, respectively) 
or of ~ 300 km s ' for the Centaurus cluster (Tamura et al. 
2014b). Only one case (Abell 2256, Tamura et al. 2011) has 
a robust detection a bulk motion of 1500 + 300 km s * been 
found. In this case, the redshift difference is consistent with 
that measured from optical spectroscopy between two sub¬ 
groups of member galaxies, and therefore the ICM is found to 
be moving together with galaxies as a substructure within the 
cluster. To summarize, very few cases of bulk motions in the 
ICM have been successfully detected so far, and most of these 
have been due to the relative motions of different massive ha¬ 
los not yet merged and still carrying their ICM, as opposed to 
bulk motions occurring within the ICM of a single virialized 
halo. 

Among the best-studied major mergers, IE 0657-56 (the 
so-called “Bullet cluster” ) occupies a special place thanks 
to the extensive investigations devoted in the last 10 years 
devoted to its dynamical structure. X-ray observations of 
ASCAand ROSAT showed a very high ICM temperature of 
~ 17 keV, and a double peaked X-ray surface brightness dis¬ 
tribution (Tucker et al. 1998). Using more accurate observa¬ 
tions from Chandra, Markevitch et al. (2002) obtained a lower 
temperature of ~ 14.8 keV. A viable explanation for such a 
high temperature is a recent, supersonic merging of two com¬ 
ponents. One of these components, called the ’’bullet”, has 
already crossed the main cluster and formed a bow shock, 
which is clearly visible in the X-ray image. The temperature 


map of IE 0657-56 is shown to be consistent with this pic¬ 
ture (Million & Allen 2009). Optical observations showed 
that the spatial distribution of the cluster’s member galax¬ 
ies was offset from the X-ray emission of the ICM (Barrena 
et al. 2002), and the gravitational potential, reconstructed with 
gravitational lensing, is also displaced with respect to the X- 
ray emission. On the other hand, the gravitational potential is 
spatially consistent with the distribution of the member galax¬ 
ies, as expected if the potential well is dominated by a colli¬ 
sionless dark matter component (Clowe et al. 2004, 2006). 
This observation has been considered as a direct confirmation 
of the existence of dark matter (Bradac et al. 2006; Angus 
et al. 2007; Clowe et al. 2007). The discontinuity in the gas 
density and temperature across the shock can be used to esti¬ 
mate the velocity of the bullet perpendicularly to the line of 
sight (Markevitch et al. 2004). The Mach number of the bow 
shock was constrained to be M = 3.0+0.4, which corresponds 
to a tangential shock velocity of ~ 4700 km s“' (Markevitch 
et al. 2006; Markevitch & Vikhlinin 2007). The radial com¬ 
ponent of the velocity was determined to be ~ 600 km s * 
based on the optical spectroscopic redshifts of member galax¬ 
ies (Barrena et al. 2002). Since this value is much lower than 
the tangential component, the direction of the bullet is ex¬ 
pected to be almost perpendicular to the line of sight. 

Although the velocity difference along the line of sight be¬ 
tween the two dark matter halos involved in the merger is 
expected to be modest (compared to the typical CCD reso¬ 
lution), we argue that the violent merger may have pushed 
a significant amount of ICM perpendicular to the direction 
of the merger, possibly resulting in significant bulk motions 
along the line of sight. In fact, we do not expect to measure 
in the ICM the same velocity value of the infalling clump, but 
instead the induced bulk motions which eventually will dis¬ 
sipate the coherent kinetic energy of the impact into thermal 
energy and turbulence in the ICM. Therefore, in this work, 
we intend to exploit the exquisite angular resolution of the 
Chandra satellite to perform spatially resolved spectral anal¬ 
ysis in several different regions to investigate the presence of 
bulk motions in the ICM of the Bullet cluster. Chandra is also 
ideal to perform this study thanks to its robust calibration and 
its well-understood background. 

The paper is organized as follows. In Section 2, we describe 
the strategy employed to search for bulk motions in the ICM. 
In Section 3, we describe a series of spectral simulations used 
to support our adopted analysis strategy. In Section 4, we de¬ 
scribe the X-ray data reduction and analysis for the entire set 
of observations of the Bullet cluster, and finally provide a dis¬ 
tribution of X-ray redshift associated with ~ 23 independent 
regions of its ICM. We will provide an accurate evaluation 
of the errors on each redshift measurement with a particular 
emphasis on the treatment of systematics associated with the 
unknown temperature structure of the ICM. In Section 5, we 
investigate the distribution of the redshift of the ICM across 
the Bullet cluster to infer the presence of bulk motions and 
discuss our results in the context of the global dynamics of 
the Bullet cluster. In Section 6, we evaluate a posteriori the ef¬ 
fect of uncertainties in the gain calibration of Chandra ACIS- 
I, through the measurement of the position of two prominent 
fluorescence lines in the background. In Section 7, we discuss 
the possible extension of this work and future perspectives for 
the detection of bulk motions in the ICM with present and fu¬ 
ture X-ray facilities. Einally, in Section 8, we summarize our 
results. Throughout this paper, we adopt the 7 year WMAP 
cosmology with 0.272, Qa = 0.728 and Hq - 70.4 km 
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s ' Mpc * (Komatsu et al. 2011). Quoted error bars always 
correspond to Icr confidence levels. 

2. A METHOD TO SEARCH EOR BULK MOTIONS IN THE ICM 

Our method to search for bulk motions in the ICM is a two- 
step process. The first step consists of identifying the pro¬ 
jected regions where the average X-ray redshift will be mea¬ 
sured. This step must be kept as simple as possible to avoid 
the introduction of bias in the redshift measurements. In fact, 
the ab initio use of spectral information, such as identifying 
regions on the basis of their hardness ratio, may lead to the 
selection of regions with maximally different spectra, which 
can, in turn, amplify the redshift differences in a non-trivial 
way. For these reasons, we will choose the independent re¬ 
gions to be analyzed only on the basis of the surface bright¬ 
ness contours, after setting an appropriate requirement on the 
minimum signal-to-noise ratio (S/N) in each region. 

The second step is to perform a standard spectral analysis to 
measure simultaneously the ICM temperature, heavy element 
abundances. X-ray redshift, and corresponding errors in the 
selected regions. This step is critical because of the unknown 
temperature structure of the ICM along the line of sight. We 
know that a significant range of temperatures is likely to be 
present in all of our spectra. We also know that the high tem¬ 
peratures involved make it impossible to constrain the differ¬ 
ent thermal components, or, in other words, the emission mea¬ 
sure of the ICM as a function of the temperature. The main 
reason is that for temperatures above 3-4 keV, any mix of dif¬ 
ferent temperatures is fit with a very good approximation by 
a single-temperature thermal emission model (Mazzotta et al. 
2004). Nevertheless, the centroid of the iron line complex de¬ 
pends significantly on the mix of temperatures involved. To 
properly take this effect into account, we choose to describe 
our spectra always as a combination of two thermal compo¬ 
nents with different temperatures. Eventually, to include the 
effects of any possible mix of different temperatures in the 
ICM, we will consider any combination of temperature pairs 
allowed by the data to evaluate the maximum variation in the 
X-ray redshift associated with the underlying temperature dis¬ 
tribution. This step allows us to fully include the effects of 
the ICM thermal structure without resorting to complex multi¬ 
temperature models, and directly provides a conservative es¬ 
timate of the systematic uncertainty on the measured X-ray 
redshift. This systematic error will be added independently to 
the statistical error. The details of this step will be given in 
Section 4.3. 

Another relevant aspect is the choice to limit our spectral 
analysis to the hard (2.0-10 keV) band. The main reason for 
this is that we prefer to base our result uniquely on the cen¬ 
troid of the iron line emission complex in the 6.7-6.9 keV rest 
frame. In principle, the inclusion of the soft 0.5-2.0 keV band 
offers the possibility to measure additional L-shell lines from 
iron and other elements, and to obtain a better measurement of 
the continuum, which is also useful to accurately measure the 
line positions. However, the inclusion of the soft band makes 
the X-ray fit much more complicated due to the possible pres¬ 
ence of small amounts of much colder gas with potentially 
different dynamical status contributing L-shell lines. There¬ 
fore, the choice of the hard band is preferred since it would 
minimize the impact of systematics associated with complex 
dynamical and thermodynamical structures. Clearly, the con¬ 
sequence is larger statistical errors with respect to the case of 
a full band analysis. As we will show with simulations, this 
effect is not large in the case of the Bullet cluster since the 


emission is dominated by very hot gas, whose emission lines 
in the soft band have a negligible contribution. In addition, 
the continuum in the soft band is strongly dependent on the 
value assumed to describe the Galactic absorption. In turn. 
Galactic absorption has a negligible impact in the hard band. 
Therefore, limiting our spectral analysis to the hard band only 
allows us to avoid another source of systematic uncertainties. 
Nevertheless, we still allow for a range of values in the Galac¬ 
tic hydrogen column density NH. 

To summarize, the spectra are fit with Xspec vl2.8.2 (Ar- 
naud 1996). To model the X-ray emission, we use two 
mekal plasma emission models (Mewe et al. 1985, 1986; 
Kaastra 1992; Liedahl et al. 1995), which include thermal 
bremsstrahlung and line emission with abundances measured 
relative to the solar values of Asplund et al. (2005) in which 
Fe/H = 3.6 X 10 ^ We note that Million & Allen (2009) 
found evidence for a power-law component in some regions 
of the Bullet cluster, possibly associated with Inverse Comp¬ 
ton on the CMB photons through a relativistic population of 
electrons. We do not attempt to model such a non-thermal 
emission, and we argue that this component, if present, would 
be accounted for by a high-temperature, low-metallicity com¬ 
ponent, and therefore its effect would automatically be in¬ 
cluded in our estimate of the total error. Galactic absorp¬ 
tion is described by the model tbabs (Wilms et al. 2000). 
The central value for the Galactic HI column density is set 
to NHcai = 4.88 X 10^°cm“^ based on Kalberla et al. (2005). 
In Million & Allen (2009), the value of NHcai at the posi¬ 
tion of the Bullet cluster has been measured from the X-ray 
analysis to be 4.44 or 4.99 x 10^° cm“^ depending on the fit¬ 
ting model, with a 1 cr error bar of 3%. We choose to allow 
NHcai to vary by an interval of 2cr (~ 6%) around the value 
of Kalberla et al. (2005). During the minimization procedure, 
the redshift parameters of both thermal components are linked 
together, while the two temperatures and the two abundances 
are left free. We apply Cash statistics (Cash 1979) to the un¬ 
binned source plus background spectra in order to exploit the 
full spectral resolution of the ACIS-I instrument. Cash statis¬ 
tics is preferred with respect to the canonical analysis of 
binned data, particularly for faint spectra (Nousek & Shue 
1989). Since it is important for us to avoid local minima 
in the fitting procedure, the fit is repeated several times be¬ 
fore and after running the command steppar on all of the 
free parameters, and particularly on the redshift, with a step 
6z = 10“^. Finally, the plots of ACstat versus redshift are visu¬ 
ally inspected to investigate whether there are other possible 
minima around the best-fit values (see Yu et al. 2011). 

3. SPECTRAL SIMULATIONS 

We run a set of spectral simulations to provide support to 
the analysis strategy outlined in Section 2, and to investigate 
our ability to recover the correct input redshift with this strat¬ 
egy. As already mentioned, the most relevant aspect we want 
to tackle when measuring the X-ray redshift is the presence 
of a multi-temperature structure along the line of sight. The 
main point is the dependence of the centroid of the iron line 
complex on the mix of temperatures present in the ICM. In 
general, we will not be able to resolve the thermal structure of 
the ICM, particularly at the high temperatures which are often 
found in massive, merging clusters. 

To investigate this aspect, we perform a first series of sim¬ 
ulations of spectra modeled with three thermal components. 
Typically, we choose two components with high tempera¬ 
tures (> 5 keV) and comparable emission measures, and a 
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Fig. 1. — Upper left panel: distribution of the best-fit X-ray redshift zx obtained fitting with a single-temperature mekal model a total of 10'* spectra simulated 
assuming a two-temperature mekal model as an input with kT[ = 8 keV and kTz = 15 keV, Zi = 0.6Zq, Z2 = 0.4Zq, and the same emission measure. A third 
component with kT^ = 3.5 keV and Z3 = 1.0 Zq has been added with a lOx lower emission measure. The redshift Zi„put = 0.30 (shown as a vertical green line) is 
the same for both components. Simulated spectra have an average of 2 X 10'* net counts in the 0.5-10 keV band. Fits ai'e perfoimed over the entire energy range 
(0.5-10 keV). Upper right panel: same as the upper left upper panel, but using a two-temperature mekal model to fit the simulated spectra. Lower left panel: 
same as the upper left panel, but fitting only the hard band (2.0-10 keV). Lower right panel: same as the upper right panel, but fitting only the hard band (2.0-10 
keV). 


colder component (~ 3.5 keV) with an emission measure 
equal to 10% of that of the hotter components. The redshift is 
Zinput = 0.30 and NHcai = 4.89 x 10^^ cm“^ for all of the com¬ 
ponents. Typical abundances are set to solar and half solar 
for the cold and hot components, respectively (in units of As- 
plund et al. 2005). The normalizations are set in order to have 
~ 20,000 net counts in the 0.5-10 keV band for an exposure of 
~ 500 ks to reproduce at best the typical situation of the actual 
data (see Section 4.1). Each set is simulated 10,000 times, and 
the background, RMF, and ARE files are taken from the real 
Chandra ACIS-I data to reproduce the instrumental set up we 
have for our spectra. As described in Section 2, in the spectra 
analysis, we link together the redshift of the thermal compo¬ 
nents while the temperatures, abundances, and normalizations 
are left free. 

We are aware that the discrete distribution of the emission 
measure as a function of the temperature in the simulated 
spectra is not meant to provide a realistic description of the 
actual ICM thermal structure; nevertheless, it is representa¬ 
tive of the situation often encountered in the ICM regions of 
the Bullet cluster, and therefore represents a useful testbed 
for our analysis strategy. In fact, it is very hard to propose a 
realistic, continuous distribution of emission measures in the 
Bullet cluster, an aspect which should be treated in its full 
complexity with hydrodynamical simulations. 

There are three relevant choices in our strategy: the method 
to find the best-fit parameters, the energy band to be con¬ 
sidered, and the fitting model. First, we set up the optimal 


TABLE 1 

Results from the simulations shown in Figure 1. 


Model 

Band (keV) 

<zx> 

<^rms 

mekal 

0.5-10 

0.3019 

0.0092 

mekal 

2.0-10 

0.3032 

0.0099 

2T mekal 

0.5-10 

0.3008 

0.0113 

2T mekal 

a 

2.0-10 

0.3013 

0.0119 


^In column 3, we show the average value of zx^ while in column 4 we show 
the dispersion cTrms of the best-ht zx distribution. 


method by running the standard analysis on our simulated 
spectra and comparing the results with those obtained with 
a simple one-step, direct fit. In our standard analysis, we re¬ 
peat the fit several times after running the steppar command 
on the redshift parameter in an interval +0.05 around the input 
redshift of Zinput = 0.3 with a step of dz - 0.0001. This proce¬ 
dure avoids the local minima and it is able to find the global 
minima in most cases. We find that in many cases (at least 
30%), the best-fit redshift obtained with a direct fit (i.e., run¬ 
ning the fit command only once) provides incorrect results 
due to a local minimum far from the input redshift. On the 
other hand, with our procedure, the correct minimum is found 
in the large majority of the cases. Therefore, we do not dis¬ 
cuss this aspect further, and henceforth adopt as the standard 
procedure the use of the steppar command with 6z - 0.0001 
on the redshift parameter. 
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TABLE 2 

Grid of temperatures and abundances used as input for the simulations 
SHOWN IN Figure 2. 


kTi 

kT2 

Z 1 /Z 0 

Z 2 /Z 0 

(o' stat) 

O' rms 

0'rms/(0' stat) 

4.0 

8.0 

1.0 

0.6 

0.0029 

0.0033 

1.13 

5.0 

lO.O 

0.9 

0.5 

0.0037 

0.0042 

1.14 

6.0 

12.0 

0.8 

0.4 

0.0048 

0.0065 

1.35 

7.0 

14.0 

0.7 

0.4 

0.0061 

0.0082 

1.35 

8.0 

16.0 

0.6 

0.35 

0.0078 

0.0108 

1.39 

#.o 

18.0 

0.5 

0.3 

0.0105 

0.0130 

1.23 


^In column five, we show the average statistical en'or on zx directly ob¬ 
tained with Xspec, while in column six we show the dispersion (Trms for tho 
distribution of zx- The ratio of these two quantities is shown in the last col¬ 
umn. 

Next, we focus on the effect of using two mekal models 
instead of one when fitting the simulated spectra, and of con¬ 
sidering only the 2.0-10 keV band instead of the entire 0.5- 
10 keV band. In Figure 1 we have an example of the ef¬ 
fect of the energy band and of the adopted model, obtained 
with input parameters kTi =8 keV, kT 2 - 15 keV, and 
Zi = 0.6Zo and Z 2 = 0.4Zo, plus a third thermal component 
with kT-j - 3.5 keV and Z 3 = I.OZq. We plot the distribution 
of the best-fit values of the X-ray redshift zx obtained while 
fitting the simulated spectra with one (left panels) and two 
(right panels) thermal components, on the 0.5-10 keV (upper 
panels) and the 2.0-10 keV (lower panels) energy bands. We 
find that in all of the cases, we are able to recover the in¬ 
put redshift with a good accuracy. However, the use of the 
two-temperature mekal model provides average redshift val¬ 
ues slightly closer to the input value, as shown in Table 1, 
where we summarize the results of this particular set of sim¬ 
ulations. Despite being very small, this effect is found in all 
of our simulations, irrespective of the input temperatures and 
abundances. We also estimate the dispersion of the zx distri¬ 
bution cTjnjs s I,i(zx,i - (zx})^/(N - 1). As expected, we find 
a larger cTj-ms when using the two-temperature model due to 
the degeneracy between the temperature values that affect the 
redshift. This effect provides ~ 20% larger errors in all of 
our simulations. On the other hand, we note that the uncer¬ 
tainties on the redshift found using the total and hard bands 
are comparable. This is due to the fact that we simulated 
high-temperature spectra, and therefore we typically have no 
relevant lines in the soft band which can add significant in¬ 
formation. We conclude that the conservative choice of using 
of a two-temperature mekal model on the hard 2.0-10 keV 
provides robust and accurate results. 

We now use simulations to investigate the accuracy of the 
statistical error provided by Xspec compared to the actual 
dispersion of the distribution of the best-fit zx- We run an¬ 
other set of simulations and analyze them directly with a two- 
temperature thermal model. Each simulation set includes 10^ 
spectra with an average of 2 x lO'* net counts as observed 
in an exposure of 500 ks with ACIS-I, with an input model 
provided by two mekal components with the same emission 
measure. The input parameters for the temperatures and abun¬ 
dances are shown in Table 2. The input redshift is always 
Zinput = 0.3. This grid has been chosen in order to represent 
the values expected in the complex temperature structure of 
the Bullet cluster (see Million & Allen 2009). For each set 
of parameters, we obtain the distribution of best-fit values of 
Zx with our standard analysis procedure, which is shown in 
Figure 2. We also derive two additional information: the av¬ 
erage 1 cr statistical error obtained with the Xspec command 


TABLE 3 

List of the 9 Chandia observations used in this work. 


ObsID 

Exposure (ks) 

Observation Date 

3184 

87.4 

2002 Jul 12 

5355 

27.4 

2004 Aug 10 

5356 

97.1 

2004 Aug 11 

5357 

79.0 

2004 Aug 14 

5358 

32.0 

2004 Aug 15 

5361 

82.6 

2004 Aug 17 

4984 

76.1 

2004 Aug 19 

4985 

27.4 

2004 Aug 23 

49k6 

41.4 

2004 Aug 25 


^All the obseiwations are taken in the VFAINT mode. In the second column 
we list the elTective exposure time after data reduction. 


“error 1 . 0 ,”, computed as (cTstat) = ^iO-smjIN, where o-stat,; 
is simply the average of the lower and upper 1 cr error bars 
on Zx estimated in the i'* simulated spectrum. Then, we also 
compute the rms dispersion of the histogram distribution as 
0 - 2 ^^, = I.iizx,i - Zinput)^/A. By comparing (cTstat) and cr^m, 
we are able to verify whether the statistical error is represen¬ 
tative of the true uncertainty on zx- In the last columns of 
Table 2 we show the values of ((Xstat), cr^jns and their ratio. 
In all of these cases, we note that is larger than the di¬ 
rect estimate of the statistical error by 15-35%. Therefore, we 
conclude that the actual uncertainty on the best-fit zx is sig¬ 
nificantly larger than the estimated statistical error. This effect 
is due to the strong degeneracy between the parameters of the 
fitting model. This is confirmed by a broad, highly scattered 
correlation between the deviation of zx with respect to the in¬ 
put value and the global deviation from the input temperatures 
that we observe in our simulations. We are also aware that the 
ratio crrn,s/(crstat) in real clusters may be even larger than the 
values listed in Table 2, since the actual thermal structure of 
the ICM is expected to be more complex than a simple two- 
temperature thermal bremsstrahlung. One last point we can 
take away from Figure 2 is that the distribution of best-fit zx 
is wider at larger temperatures. This simply reflects the dif¬ 
ficulty in measuring the iron emission-line complex in high- 
temperature spectra. 

Taken together, these results support our choice of mea¬ 
suring an additional systematic error associated with the un¬ 
known thermal structure of the ICM, as described in Section 
2. To do that, for each region, we will derive the best-fit X-ray 
redshift for all of the possible combinations of temperatures 
and abundances sampled in a realistic range (from 3.5 to 27 
keV and from 0.1 to 1.6 Zq, respectively) and measure the 
width of the distribution of the best-fit Zx(Ti, 7 ’ 2 ,Zi,Z 2 ). This 
will provide us with the systematic uncertainty cTsyst on zx 
associated to the unknown thermal composition of the ICM. 
The details of this step and the selection of a proper range of 
temperatures and abundances in each region will be described 
in full detail in Section 4.3. Finally, this uncertainty will be 
added to the statistical error to provide the total uncertainty 
CTtot on the best-fit zx- 

4. X-RAY DATA REDUCTION AND ANALYSIS 
4.1. Data reduction 

The Chandra observations used in this paper are listed in 
Table 3. We find 10 pointings of IE 0657-56 in the Chan¬ 
dra archive, of which we choose 9 after discarding Obsid 554 
which was taken in the FAINT mode. Due to the low expo¬ 
sure of Obsid 554 (~ 26 ks) this choice has little effect on 
our final results. The selected observations were carried out 
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Fig. 2. — Histogram distribution of the best-fit redshift zx for 10^ simulated spectra with input temperature and abundances listed in Table 2 and Zinpm = 0.3 
(shown with a green vertical line). In red, we show the Gaussian with amplitude crpms centered on (zx)- 


between 2002 July and 2004 August in VFAINT mode using 
the ACIS-I. The data reduction is performed using the latest 
release of the ciao software (version 4.6) with CALDB 4.6.7. 
Charge transfer inefficiency (CTI) correction, time-dependent 
gain adjustment, grade correction, and pixel randomization 
are applied. We are able to filter efficiently the background 
events thanks to the VFAINT mode. This allows us to reduce 
the background by a ~ 25% in the hard (2.0-10 keV) band. 
Eventually, we search for high background intervals and re¬ 
move them with a 3-cr clipping. The final exposure times of 
each Obsid, listed in Table 3, are lower than the nominal ex¬ 
posure time only by a few percent. The level 2 event files 
obtained in this way are then reprojected to match the coordi¬ 
nates of Obsid 3184, and merged to a single event file. The 
total exposure time of the merged data is ~ 550.4 ks. The 
soft- (0.5-2.0keV), hard- (2.0-lOkeV) and total- (0.5-10keV) 
band images are obtained from the merged file. 

We select the regions from which we extract the spectra on 
the 0.5-10 keV band image. We apply the contour-binning 
technique of Sanders (2006), developed in order to select re¬ 
gions according to the surface brightness distribution of an 
extended source. Our goal is to obtain spectra with compara¬ 
ble constraints on the X-ray redshift. To do that, we ideally 
should require that there be a similar number of net counts in 
the 2.0-10 keV band in each region. However, since the range 
in region size is wide and we have a few regions with a much 
larger area than average (particularly in the outer part of the 
clusters), we prefer to set our criterion in terms of S/N in the 
total band. Therefore, in this way we select 23 regions with 
full-band S/N> 140 each within a radius of 104", which is the 


circular region with the highest S/N in the 0.5-10 keV band, 
and will be used to compute the global redshift of the ICM 
(see Section 4.2). 

In Table 4, we list the net counts measured in the soft and 
hard bands, and the S/N in the hard band, for the 23 regions 
defined by the contour-binning algorithm. We have at least ~ 
7000 net counts in each region in the 2.0-10 keV band with a 
maximum of 12,000/14,000 net counts for the largest regions. 
Note that this value is much larger than that generally used 
to measure temperature, iron abundance, and average X-ray 
redshift. Here, we set a higher threshold because we want to 
achieve a precision on the X-ray redshift measurement larger 
than usually needed to evaluate global properties, and because 
of the exceptionally high temperatures in the ICM that make 
the iron emission-line complex less prominent. The discrete, 
color-coded surface brightness of the Bullet cluster, together 
with the shape of the selected regions, is shown in Figure 3. 

Finally, we extract the spectra of each region in the full 
band. Thanks to the exquisite angular resolution of Chan¬ 
dra, we do not need to correct for effects due to the point 
spread function when extracting the spectra. Before creat¬ 
ing the spectra, we manually remove all of the point sources 
visible in the soft- and total- band images. For each region, 
we produce response matrix (RMF) and ancillary response 
(ARF) files from each Obsid, and then add the RMF and ARF 
files, weighing them by the corresponding exposure times. In 
this way, we keep track of all the differences in the ACIS- 
I effective area among different regions on the detector and 
among Obsid taken at different epochs. This is an important 
step, considering that the cluster is observed on different CCD 
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Fig. 3.— Color-coded surface brightness map of the Bullet clusters. Each of the 23 regions defined by the contour-binning technique are shown with a different 
color corresponding to the hard-band surface brightness. The units of the color bar are counts per pixel. The large circle shows the region used to extract the 
spectrum of the global cluster emission. 


TABLE 4 

Soft- and hard- band photometry of the 23 regions selected for spectral 

ANALYSIS. 


Reg ID 

Net Cts (0.5-2.0keV) 

Net Cts (2.0-lOkeV) 

S/N (2.0-lOkeV) 

0 

13402 ± 116 

6519 ±82 

78.8 

1 

12322 ± 111 

7484 ± 88 

84.1 

2 

12678 ± 113 

7049 ± 86 

81.1 

3 

12355 ± 112 

7410± 88 

83.4 

4 

20335 ± 149 

12432 ± 133 

93.3 

5 

12306 ± 111 

7418 ±89 

83.1 

6 

12267 ± 111 

7412 ± 89 

82.8 

7 

12316± 112 

7185 ± 88 

80.8 

8 

12138± 111 

7384 ± 89 

82.1 

9 

12009 ± no 

7627 ± 91 

83.5 

10 

12176± 111 

7266 ± 90 

80.6 

11 

15556 ± 129 

9670 ± 113 

85.2 

12 

12595 ± 115 

7946 ± 99 

79.7 

13 

12274 ± 112 

7397 ± 90 

81.6 

14 

11928 ± 111 

7374 ± 91 

80.6 

15 

11946± 110 

7586 ±91 

82.8 

16 

11900± 111 

7142 ±91 

78.0 

17 

20063 ± 144 

12160 ± 120 

100.8 

18 

21797± 153 

14097 ± 136 

103.4 

19 

11691 ± no 

7496 ± 93 

80.0 

20 

11880 ± 111 

7250 ± 92 

78.1 

21 

12068 ± 111 

7214 ± 90 

79.5 

22 

11632 ± 109 

7665 ± 92 

82.5 


^Errors on the net detected counts are computed as Vo^ + 2x5, where 
Cts are the source counts and B is the expected background geometrically 
rescaled to the source region. 


(ccd_id=3 for Obsid 3184, 4984, and 5361; ccd_id=2 for 
Obsid 4985,5355, and 5356; ccd_id=® for Obsid 4986,5357, 
and 5358). The background is sampled from a series of cir¬ 
cular regions as far as possible from the cluster emission, but 
still on the solid angle dehned by the overlap of all the Obsid. 
In this way, we are not sampling regions where the exposure 
time is lower than 550 ks, and therefore the background can 
be simply geometrically scaled to the source region. The geo¬ 
metrical backscale values between the background region and 
the source region range from ~ 10 to ~ 100, and the ratio of 
the background to the ICM signal in the source region ranges 
from 2% up to 15% for the outermost region. 

4.2. Global ICM redshift 

As a hrst step, we measure the redshift of the ICM emis¬ 
sion from the entire cluster, htting the global spectrum in the 
2.0-10 keV energy range. The extraction region of the global 
ICM emission is dehned as the circle with the maximum S/N 
(see Yu et al. 2011) corresponding to a radius of 104”. This 
circular region also provides the outer limit for the dehnition 
of the ICM regions, as shown in Figure 3. 

Given the complex temperature structure encompassed by 
this large region, we hnd that one of the two temperatures 
hits the hard upper limit during the minimization proce¬ 
dure. Therefore, we freeze it to a reference maximum value 
of 26 keV. The best-ht X-ray redshifts obtained from the 
spectral analysis of the global emission from the cluster is 
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Zxgiobai = obtained for two temperatures of 

26 keV (frozen) and 6.2;|;[ 3 keV, and iron abundances of 
(0.50 + 0.10 )Zq and OAQ^qq^Zq for the first and the second 
thermal components, respectively. We note that the best-fit 
value of the Zxgiobai is consistent within 1 cr with the optical 
value Zopt = 0.296. We remark that we would obtain a value 
of 0.3059 + 0.005, therefore, at ~ 2cr from the optical value 
if we had used a single mekal model instead of two. This is 
another hint toward the necessity of using at least two tem¬ 
peratures to properly model the thermal structure of the ICM, 
and therefore the shape of the iron line complex. We also note 
that the best-fit Zxgiobai is inconsistent with the value z = 0.325 
published in Yu et al. (2011). This is mostly due to the fact 
that in Yu et al. (2011), we used a one-temperature mekal 
model and adopted a much simpler minimization procedure. 
Our aim was to test the capability of recovering the X-ray red- 
shift with a blind spectral analysis with a single-temperature 
mekal model, since this procedure was meant to be applied 
to a poorly characterized cluster. Therefore, the discrepancy 
found in Yu et al. (2011) is affected by the temperature struc¬ 
ture of the ICM. 

This result shows that the global ICM emission is dom¬ 
inated by gas which is moving along with those galaxies, 
whose velocities lie mostly along the plane of the sky, with 
a small velocity component projected along the line of sight 
estimated to be ~ 600 km s * based on optical spectroscopy 
(Barrena et al. 2002). This velocity component would cor¬ 
respond to a maximum difference between the galaxies and 
the majority of the ICM of Az ~ 2.6 x 10“^, which is be¬ 
low the detection threshold of our analysis. As expected, the 
only means to detect any internal bulk motions induced by the 
merger is to restrict our analysis to small regions, exploiting 
the exquisite angular resolution of Chandra to search for re¬ 
gions where a relevant mass of gas has reached a velocity of at 
least few x 1000 km with respect to the bulk of the ICM. 

4.3. Spatially resolved spectral analysis with contour binning 

In this Section, we present the results from the fits of the 
23 independent regions selected by the contour-binning tech¬ 
nique. Fits are performed in the 2.0-10 keV energy range ac¬ 
cording to the strategy provided in Section 2. For each region, 
we produce a plot showing the ACstat value versus the redshift, 
obtained by varying the redshift parameter and marginalizing 
the fit with respect to the other parameters. Two examples are 
shown in Figure 4. The minimum is well defined and roughly 
symmetric, therefore providing a robust estimate of the sta¬ 
tistical uncertainty. Nevertheless, we visually inspected all of 
the ACstat versus redshift plots, finding one anomaly in the 
case of regions 7 and 15, where two comparable minima are 
found. This may be due to the overlap of two comparable 
masses of ICM with different redshift, but we cannot confirm 
this interpretation. We do not find similar cases in the remain¬ 
ing 21 regions. Therefore, we will exclude regions 7 and 15 
when evaluating the probability of ICM bulk motions in Sec¬ 
tion 5. 

In Table 5, we list the best-fit zx with the corresponding 
lower and upper statistical error bars corresponding to the Icr 
confidence level. This list of values defines a redshift map of 
the ICM in the Bullet cluster, shown in Figure 5. The map 
appears to be patchy with the strongest redshift differences 
in the trail of the bullet, while there is no striking difference 
between the bullet and the surrounding gas, consistently with 
the finding that the merger is occurring in the plane of the sky. 




Fig. 4.— Plot of ACstat vs. redshift for region 4 (left panel) and region 9 
(right panel). 

The Zrmx map suggests that significant bulk motions may be 
left behind the bullet’s path, which eventually will evolve into 
turbulent motions. 

The redshift map has also been derived using another back¬ 
ground, chosen from an independent (not overlapping) region 
on the detector with the same characteristics as described in 
Section 4.1. We find minor differences which leave our re¬ 
sults unaltered. We also repeated the fit with a synthetic back¬ 
ground obtained from blank field Chandra data after normaliz¬ 
ing it in the 2.0-10 keV band to our reference background ex¬ 
tracted from the data. Again, we find consistency well within 
1 cr for the best-fit parameters. Therefore, we conclude that 
our results are robust against variation in the background. 

Before drawing any conclusion on any possible bulk mo¬ 
tion of the ICM caused by the passage of the bullet, we need 
to carefully take into account the uncertainty on the velocity 
measurement in each region. However, as discussed in Sec¬ 
tion 2 and shown with simulations in Section 3, the statistical 
errors may underestimate the actual uncertainties on the best- 
fit Zx, and therefore on the relative velocity in the observed 
frame vicm - cx {zi - (zx))/(l + {zx})- Therefore, we need 
to evaluate in each region the uncertainty associated with the 
unknown thermal structure of the ICM along the line of sight. 

To do so, we proceed as described in Section 3. We con¬ 
sider a four-dimensional grid of spectral parameters Ti, T 2 , 
Z\ and Z 2 covering all of the possible values expected in the 
Bullet cluster. We set a lower limit of 3.5 keV and a maximum 
of 27 keV for the temperatures. As for the metal abundances, 
we make them range from 0.1 to 1.6 times the solar value 
(in units of Asplund et al. 2005). The steps of the grid are 
0.5 keV for the temperature and 0.1 for the metal abundance. 
Clearly, this temperature and abundance grid is exceedingly 
wide for a single ICM region. Therefore, we selected a sub¬ 
grid in the following way. First, for each region, we measure 
the best fit in the full (0.5-10 keV) band, leaving all of the pa¬ 
rameters free, and collect the absolute minimum Cmin. Then, 
we compute the best fit on the grid in the full (0.5-10 keV) 
band where only the redshift parameter and the two mekal 
normalizations are left free, while NHcai is allowed to vary 
in a 6% interval around the central value. Thus, we com¬ 
pute the best fit achievable for each set of values on the grid, 
and finally select all of the values whose best-fit Cstat is close 
enough to the absolute minimum Cmin- The criterion we adopt 
is given by ACstat = Qtat - Cmin < 4.72, which corresponds 
to a Icr confidence level for four free parameters (see Press 
et al. 2002). In this way, we also exploit also the information 
included in the soft band (0.5-2 keV). By applying the crite¬ 
rion ACstat < 4.72 for the fits performed in the full 0.5-10 keV 
band, we select the subgrid of Ti, 7’2, Zi, and Z 2 values statis¬ 
tically compatible within 1 cr with the spectrum observed in 
each region. Then, we again run the spectral fit on the subgrid 
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TABLE 5 

Results of the spectral fits of the ICM projected regions shown in Figure 3. 


Region 

Z 

O ' stati , 

Statu 

O ' systi , 

O'systu 

O'totb 

i^totu 

0 

0.3000 

- 0.0053 

0.0047 

- 0.0025 

0.0001 

- 0.0058 

0.0047 

1 

0.3265 

- 0.0099 

0.0090 

- 0.0050 

0.0001 

- 0.0110 

0.0090 

2 

0.2994 

- 0.0057 

0.0027 

- 0.0001 

0.0080 

- 0.0057 

0.0084 

3 

0.2957 

- 0.0159 

0.0201 

- 0.0001 

0.0169 

- 0.0159 

0.0262 

4 

0.2787 

- 0.0055 

0.0081 

- 0.0001 

0.0048 

- 0.0055 

0.0094 

5 

0.2874 

- 0.0105 

0.0087 

- 0.0029 

0.0024 

- 0.0108 

0.0090 

6 

0.2957 

- 0.0113 

0.0105 

- 0.0047 

0.0002 

- 0.0122 

0.0105 

7 

0.2836 

- 0.0054 

0.0088 

- 0.0001 

0.0072 

- 0.0054 

0.0113 

8 

0.2719 

- 0.0097 

0.0067 

- 0.0032 

0.0070 

- 0.0102 

0.0096 

9 

0.3061 

- 0.0077 

0.0108 

- 0.0026 

0.0018 

- 0.0081 

0.0109 

10 

0.3159 

- 0.0102 

0.0092 

- 0.0056 

0.0001 

- 0.0116 

0.0092 

11 

0.3062 

- 0.0132 

0.0122 

- 0.0027 

0.0007 

- 0.0134 

0.0122 

12 

0.3215 

- 0.0106 

0.0095 

- 0.0087 

0.0001 

- 0.0137 

0.0095 

13 

0.3181 

- 0.0167 

0.0137 

- 0.0099 

0.0050 

- 0.0194 

0.0145 

14 

0.2934 

- 0.0134 

0.0302 

- 0.0026 

0.0046 

- 0.0136 

0.0305 

15 

0.3083 

- 0.0041 

0.0085 

- 0.0038 

0.0001 

- 0.0055 

0.0085 

16 

0.2933 

- 0.0154 

0.0187 

- 0.0023 

0.0020 

- 0.0155 

0.0188 

17 

0.3138 

- 0.0317 

0.0181 

- 0.0179 

0.0001 

- 0.0364 

0.0181 

18 

0.3070 

- 0.0086 

0.0084 

- 0.0056 

0.0001 

- 0.0102 

0.0084 

19 

0.2665 

- 0.0241 

0.0175 

- 0.0071 

0.0034 

- 0.0251 

0.0178 

20 

0.2922 

- 0.0090 

0.0134 

- 0.0001 

0.0103 

- 0.0090 

0.0169 

21 

0.2935 

- 0.0091 

0.0104 

- 0.0001 

0.0040 

- 0.0091 

0.0111 

22 

a 

0.2941 

- 0.0195 

0.0191 

- 0.0006 

0.0026 

- 0.0195 

0.0192 


“Column 1: region number; Column 2: best-fit redshift zx obtained fitting the 2.0-10 keV energy range; Columns 3 and 4; lower and upper Icr error 
bars from fit statistics; Columns 5 and 6: lower and upper Icr error bars from systematics associated with the ICM temperature structure; Columns 7 


and 8: total lower and upper Icr error bars computed as crtot = 
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Fig. 5.— Redshift map of the ICM in the Bullet cluster computed with respect to the average redshift (zx) = 0.2986. The color bar used discrete steps equal to 
Az = 0.007 which corresponds approximately to the typical 1 cr error on zx- 
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but in the 2.0-10 keV band only, and measure a set of best- 
fit zx values corresponding to the subgrid parameters. The 
distribution of zx obtained in this way is used to derive a red- 
shift range defined by the upper and lower 90% percentiles 
of the Zx distribution on the grid. Two examples are shown 
in Figure 6, where the upper and lower 90% percentiles are 
marked with vertical dashed lines. We remark that the dis¬ 
crete appearance of the distribution is due to the spacing of 
the parameter grid, since in several cases the best-fit redshift 
is not affected by the metal abundance value. A finer and more 
time-consuming grid would provide a smoother distribution 
with little effect, however, on the upper and lower 90% per¬ 
centiles, as we verified. This range provides an estimate of the 
systematic uncertainty intervals due to the unknown thermal 
structure. We remark that the distribution of ZxiTi, T 2 ,Z\,Z 2 ) 
is already broadened by the statistical uncertainty. Therefore, 
we are confident that the lower and upper 90% percentiles are 
a conservative upper limit to the real 1 cr uncertainty due to 
the thermal structure of the ICM. 

We argue that this systematic uncertainty cr^y^i is uncor¬ 
related with the statistical error since it is due to the de¬ 
generacy of the temperatures values, as discussed in Sec¬ 
tion 3. Therefore, we estimate the total Icr uncertainty as 

f^tot “ + The values of cTsyst and (Xtot are shown 

in Table 5 where we list separately the lower and upper er¬ 
ror bars. Although the errors are not symmetric, in Figures 
7 and 8, only for the sake of visualization we show a map of 
the statistical and systematic cr error bars across the cluster 
image, averaging the lower and upper error bars. The color 
scale shows the accuracy on zx, with brighter regions having 
smaller errors. It is possible to see that the systematic uncer¬ 
tainties are typically smaller than the statistical ones. 

A possible way to combine the information in Figure 5, 7, 
and 8 is to visualize the significance of the redshift difference 
as (zx - <2x))/o'tot, which is shown in Figure 9. Here, <zx) 
is the average zx as computed in Section 5. Bright blue and 
bright red regions point out ~ 3cr deviations with respect to 
the average redshift. It is possible to note that most of the re¬ 
gions are within 2cr from {zx), while two symmetric regions 
show a 3 cr deviation, approximately around the central region 
of the Bullet cluster. In Figure 10, we show the binned spec¬ 
tra of the two regions with the most extreme redshift (region 1 
withzx = 0-3265 ^qqjjq and region 8 withzx = 0-2719 qq°q|®). 
We focus on the 4.0-7.0 keV energy band to emphasize the 
difference in the position of the iron line complex in the two 
regions. In particular, in the right panel of Figure 10, we 
show the data and best-fit model of region 8 (red solid line 
and points) along with the best-fit model of region 1 (black 
solid line). Despite the heavy binning with S/N=9 per bin, 
it is possible to visually appreciate the difference in the posi¬ 
tion of the iron line complex among the two regions. In the 
following section, we will further investigate and quantify the 
statistical significance for the presence of bulk motions across 
the ICM. 

5. RESULTS ON THE PRESENCE OF BULK MOTIONS 

In Figure 11, we show the histogram distribution of the 
best-fit redshift. The distribution clearly appears as non- 
Gaussian. However, it is crucial to consider the information 
on the total error before drawing any conclusions. A direct 
evaluation of the presence of bulk motion in the ICM of the 
Bullet cluster is obtained with a simple test of the distribu¬ 
tion of Zx- We consider all of the regions, except for regions 


7 and 15, and Und x^ = 2(zx - {zx))^/o'^ot - 32.2, which, 
for 20 degrees of freedom, corresponds to a one-tailed (right- 
tail) probability value for a x^ test of ~ 0.04. If we consider 
the value obtained with our alternative background choice, we 
findx^ = 32.3. Therefore, we can conclude that the statisti¬ 
cal evidence for inhomogeneities in the overall distribution of 
Zx among the ICM regions of the bullet cluster is significant 
at a ~ 2 cr confidence level. Currently, the main limitation 
we encounter in measuring the bulk motions in the ICM of 
the Bullet cluster arises due to the statistical error. A better 
spectral resolution would help to bring this error further be¬ 
low the systematic one, possibly allowing us to confirm the 
detection of bulk motion obtained in this work. On the other 
hand, the systematic error can be reduced only by better con¬ 
straining the high ICM temperatures along the line of sight, 
which would require a high efficiency, high-resolution, hard 
X-ray telescope, which is not currently available or planned. 

We also search for any possible correspondence between 
the ICM redshift map and the redshift distribution of the 
member galaxies. In the SIMBAD astronomical database, we 
search for galaxies with spectroscopic redshift and within 20 
arcmin of the cluster center. We find 108 galaxies meeting 
these conditions. Then, we fit the distribution of velocities 
with a Gaussian and perform a 3cr clipping by iteratively re¬ 
moving the galaxies whose velocity lies outside +3cr from 
the mean velocity. We finally identify 80 member galaxies, 
in agreement with previous works by Barrena et al. (2002) 
and Guzzo et al. (2009). We compute the mean velocity of 
the member galaxies as Vd = 88766 + 1252 km s“' based 
on the Gaussian best fit, corresponding to an optical redshift 
of zo=0.296+0.004. Both values are consistent with Barrena 
et al. (2002) and Williamson et al. (2011). We then compute 
the peculiar velocities of each member galaxy, finding a range 
from -3300 km s"' to 3200 km s“^ 

We overplot the member galaxies and their rest-frame ve¬ 
locities on the redshift map (see Figure 12). Each circle with a 
dot or cross in its center represents a galaxy; the radius of the 
circle is proportional to the peculiar velocity of the galaxy. 
Red circles and red crosses stand for member galaxies with 
positive velocity along the line of sight, while member galax¬ 
ies with negative velocities are represented by blue circles and 
blue dots. We do not find any hint of a correlation between the 
motion of the galaxies and the velocity map of the ICM, sug¬ 
gesting that in any case there are no halos whose galaxies and 
ICM move together along the line of sight. We conclude that 
if there is a significant bulk motion in the ICM, then this is de¬ 
coupled from the galaxy motions, and therefore it is likely to 
be caused by disordered velocity fields within the ICM gen¬ 
erated by the merger event pushing away significant amounts 
of ICM mass. 

To further investigate the direction along the trail of the 
“bullet”, we select a direction marked in green in Figure 12. 
In Figure 13, we show the redshift measured in each region 
crossed along this direction. We note that approximately at 
the center of the cluster, within the box drawn in Figure 13, 
the projected redshift of the ICM appears to change signifi¬ 
cantly within 200-300 kpc, while outside this box the values 
of Zx are consistent with the average redshift, shown with a 
horizontal line. Nominally, the maximum different redshift is 
Az = 0.054 + 0.015. We can express this results by stating 
that for a region of ~ 200 - 300 kpc, we observe a velocity 
gradient of 46 ± 13 km s^'/kpc. We argue that this may be the 
signature of a significant amount of ICM pushed along the 
line of sight (perpendicularly to the merging) at a velocity of 
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Fig. 6. — Histogram distribution of best-fit zx(T\,T 2 ,Z\,Z 2 ) computed on the subgrid of temperature and abundance values selected as described in the text, for 
region 5 (left panel) and region 8 (right panel). The central vertical line shows the absolute best-fit zx^ the dashed area the range due to the statistical errors on 
the absolute best-fit zx and the two vertical dashed lines the lower and upper 90% percentiles of the zx{T[,T 2 ,Zi,Z 2 ) distribution. 
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Fig. 7.— Map of the statistical Icr error on the X-ray redshift zx- 
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Fig. 8.— Map of the systematic Icr error on the X-ray redshift zx- 


several thousand km s ' in both directions. 

To summarize, we find marginal evidence of bulk motions 
in the ICM of the Bullet cluster at a 2cr confidence level. If 
we focus on the bullet path, then we find that the maximum 
redshift difference is significant at a 3cr level between two spe¬ 
cific, symmetric regions. If confirmed, then this finding would 
imply the presence of shocks along the line of sight. We argue 
that the combination of high-resolution spectroscopy, and spa¬ 
tially resolved, moderate resolution spectroscopy may help to 
resolve the dynamical structure of the ICM in the Bullet clus¬ 
ter. 

6. EFFECTS OF GAIN CALIBRATION UNCERTAINTIES 

So far, we have ignored the effects of the uncertainties in 
the ACIS-I energy scale calibration (i.e. gain calibration). For 
example, it is known that there are shifts in the energy of the 
fluorescent emission lines in the ACIS background due to the 
change of the detector temperature (see Sanders et al. 2014). 
To investigate whether this happens in the data of the Bullet 
cluster, we accurately fit the centroids of the Ni Ka and Au La 
fluorescence lines, which are prominent in the ACIS-I back¬ 
ground spectrum at ~ 7.5 and ~ 9.7 keV, respectively. We 
simply fit in a narrow energy range (~ half a keV) a power- 
law continuum plus a Gaussian line at the position of the two 
lines. We consider all of the emission from the CCD after ex¬ 
cluding a circle of 104” centered on the cluster and all of the 
point sources. This does not remove all of the cluster emis¬ 
sion, but the fluorescence lines strongly dominate the contin¬ 
uum components associated with the ICM emission and the 
background. The best-fit energies of the two lines in each 
Obsid are shown in Figure 14. If we consider the statistical 
error associated with the line energies, then we find that the 
intrinsic shift in energies between the Obsid is ~ 0.010 keV, 
which corresponds roughly to Az ~ 0.0014. This is at least 
five times smaller than the total uncertainties in zx- Therefore, 
we conclude that the global shift in gain calibration between 
the nine Obsid does not significantly affect our results and, in 
any case, should have a mild effect in broadening the iron line 
profile observed in the total spectrum of each region. 

However, we cannot exclude the spatial variation of the gain 
calibration, which can amplify the fluctuations in the best-fit 
Zx values. To estimate the possible bias due to this effect, 
we again fit the Ni Ka and Au La fluorescence lines in the 
spectra of the 23 regions. In this case, the statistical errors 
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Fig. 9. — Significance map obtained as (zx ” Zavg)/(Tio{- Only two symmetric regions stand out in bright red and blue colors, while all of the remaining regions 
have zx within less than 2(T from (zx)- 



Fig. 10.— Folded spectra of the two regions with the most extreme redshift (left panel: region 1; right panel: region 8). The spectra are binned to a minimum 
of S/N=9, and only the 4.0-7.0 keV energy range is shown. The solid lines show the best-fit model obtained fitting the 2.0-10 keV energy band with two mekal 
components. The dashed line in the right panel shows the best-fit model of region 1. 


associated with the line centroids are much larger, due to the 
smaller statistics and the much higher contribution from the 
ICM emission above 7 keV. However, we are able to mea¬ 
sure the line centroid for Ni Ka in most of the regions and for 
Au La in all of the regions (see Figure 15). The scatter mea¬ 
sured among the regions is consistent with statistical noise. 
A simple test assuming a constant value provides a one- 
tailed (right-tail) probability value of 0.16 for both lines. This 
means that there is no evidence for significant gain-calibration 


variations from region to region that could amplify the zx fluc¬ 
tuations found in our analysis. More important, we find that 
the energy shift we measured, taken at face value (i.e., assum¬ 
ing that the energy shift is entirely due to an intrinsic change 
in the gain calibration) typically anticorrelates with the red- 
shift difference we found in our analysis, in particular, in the 
regions with the largest Az (regions 8 and 1). If we simply 
subtract the gain shift map from our zx map, then the statis¬ 
tical significance of a deviation from a Gaussian distribution 
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Fig. 11. — Histogram distribution of the best-fit reshift zx of the 23 ICM regions. The average redshift (zx) = 0.2986 is shown as a vertical line. 


increases. Clearly, we do not consider this correction in our 
final results, since the gain-fluctuation map is likely to be en¬ 
tirely dominated by the statistical noise. 

To summarize, we do not attempt to provide a full treatment 
of the gain-calibration uncertainties in ACIS-I data, but we 
perform a direct empirical test which allows us to conclude 
that any effect associated with gain-calibration uncertainty is 
safely below the typical statistical errors on our zx values, and 
therefore does not significantly affect our conclusions. 

7. DISCUSSION 

In this work, we developed a method to search for bulk mo¬ 
tions in the ICM with Chandra, the X-ray facility which is 
currently providing the best performance in terms of spatially 
resolved X-ray spectroscopy. Despite this, we are not able 
to provide a definitive evidence for bulk motions along the 
line of sight in the massive merger cluster IE 0657-56. We 
are currently preparing a systematic search of bulk motions 
in a sample of merging or post-merger clusters observed with 
Chandra, in order to achieve a stronger statistical detection of 
bulk motions. Currently, we acknowledge that this kind of 
study is challenging and may provide positive results only in 
a few extreme cases. 

Without any doubt, mapping the velocity field in the ICM 
in the majority of clusters would be a key tool to studying the 
process of cluster formation, and to achieving more accurate 
mass measurements for cosmological studies. Indeed, this is 
one of the science goals of future planned X-ray missions with 
micro calorimeters, like the SXS instrument (Mitsuda et al. 
2010) on board ASTRO-H (Takahashi et al. 2012), which 
will also reach the 6.7-6.9 keV energy range corresponding 


to the rest-frame Kg, iron emission lines. With an average 
spectral resolution of about 7 eV, it is possible to find bulk 
motions at least in the brightest nearby clusters (see Tamura 
et al. 2014a). In addition, the line broadening can be mea¬ 
sured opening the way to the investigation of ICM turbulence. 
However, a detailed study requires sophisticated approaches 
to properly measure the profiles of the emission lines (see, 
e.g., Shang & Oh 2012) to overcome the limitations due to 
the poor angular resolution. 

In the future, bulk motion studies will be possible both with 
CCD and microcalorimeter spectra thanks to the Athena mis¬ 
sion (Barcons et al. 2012). The much higher effective area 
with respect to Chandra will help to increase significantly the 
accuracy of the centroid of the iron emission lines. However, 
the limited spectral resolution of the CCD data will not help 
to remove the spectral degeneracy among different spectral 
components, and therefore the angular resolution still remains 
crucial to analyzing separately regions with different dynami¬ 
cal properties. On the other hand, the X-ray Microcalorimeter 
Spectrometer (XMS) on board Athena will be able to achieve 
a redshift accuracy of 10“^ in one pointing with an exposure of 
only 100 ks for a cluster like A2256 (see Nevalainen 2013), 
which is three orders of magnitude larger than the accuracy 
achieved in this work. Clearly, the angular resolution will be 
crucial in order to investigate the 3D dynamical state of clus¬ 
ters, particularly at high redshift where the fraction of major 
mergers is expected to increase significantly. 

X-ray spectroscopy is not the only means to investigate bulk 
motions in the ICM. In particular, the bulk motion in IE 0657- 
56 may also be detected by the kinetic Sunyaev Zeldovich ef¬ 
fect (kSZE) caused by the Doppler effect of the cluster radial 
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Fig. 12.— Spatial distribution of 80 member galaxies over-plotted on the ICM redshift map. The radius of the circles is proportional to the peculiar velocity of 
the galaxy. Red circles and red crosses stand for member galaxies with Vpest > 0, while blue circles and blue dots stand for member galaxies with Vpest < 0. The 
green arrow shows the directions we consider to explore the X-ray redshift along the bullet trail (see Figure 13). 



d(kpc) 


Fig. 13.— Redshift measured along the direction of the bullet where the 
x-axis is the distance from the peak of the X-ray brightness (the head of the 
bullet). The average redshift (zx) = 0.2986 is shown as a horizontal line. 
The dashed box roughly identifies the region in which we observe a velocity 
gradient of (46 ± 13) km s“*/kpc. 

bulk motion (see Carlstrom et al. 2002). The advent of high- 
resolution SZE imaging of clusters, reaching the 10” - 30” 
regime, opens a new window on the study of ICM bulk mo¬ 
tions. In a few cases, SZE enhancements have been inter¬ 
preted as effects of the ICM dynamics (see Mason et al. 2010, 
for the case of RXJ1347). Recent simulations have investi¬ 
gated this aspect, showing that kSZE observations may be ef¬ 
fective for probing the ICM motions in major mergers. How¬ 
ever, based on our results, the kSZE in the Bullet cluster is 


expected to appear on scales of ~ 4” - 5", which requires 
an angular resolution of at least 2", which is out of reach for 
present-day SZ telescopes (Kelley et al. 2010; Barcons et al. 
2012). Therefore, the detection of kSZE in the Bullet clus¬ 
ter appears to be challenging despite the strong SZE signal 
observed using the APEX-SZ instrument at 150 GHz with a 
resolution of 1' (Halverson et al. 2009). To summarize, the 
investigation of ICM bulk motions across the majority of the 
known cluster population will be extremely difficult in the 
near future, and a major breakthrough in this field may be 
achieved only with a specially designed instrument. 

8. CONCLUSIONS 

We propose a simple method to search for ICM bulk mo¬ 
tions in Chandra CCD data, exploiting the high angular reso¬ 
lution to perform spatially resolved X-ray spectroscopy. The 
method consists of fitting spectra extracted from the ICM re¬ 
gion selected on the basis of the surface brightness contours. 
The fitting model is given by two mekal components with 
all of the parameters set free except the redshift, which is the 
same for the two thermal components. The fits are performed 
on the hard 2.0-10 keV band only to focus on the 6.7-6.9 keV 
He-like and H-like iron emission lines. 

We apply this method to the Bullet cluster IE 0657-56, a 
bright, massive merging cluster at z ~ 0.3. The choice of this 
cluster for testing our method is based on an expectation of 
high bulk motions due to the large kinetic energy associated 
with the merger of two massive clusters. On the other hand, 
the large temperatures involved, with kT > 10 keV, make the 
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Best-fit energies for the fluorescent lines of Au La (upper panel) and Ni Ka (lower panel) lines in each Obsid with the corresponding Icr en'or bars. 


Fig. 14.— 

spectral diagnostics quite difficult for an X-ray instrument like 
Chandra ACIS-I, which has an effective area rapidly decreas¬ 
ing above 7 keV. In order to reach the required accuracy in 
redshift, we select regions with about 2x10"^ photons in the 
full 0.5-10 keV band, a requirement much stronger than that 
usually used to measure the X-ray redshift. 

We measure zx in 23 independent regions, and in each re¬ 
gion we accurately evaluate the statistical error on zx and the 
systematic uncertainty associated with the unknown thermal 
structure of the ICM along the line of sight. This allows us to 
estimate the total Icr error on zx- We find that the distribution 
of Zx is inconsistent with a constant redshift across the ICM 
of the Bullet cluster at the 2cr confidence level. If we focus on 
a particular direction corresponding to the path of the bullet, 
then we observe a velocity gradient of 46 + 13km kpc^^ 
for a region 200-300 kpc in size in the center of the cluster. 
Finally, we verify a posteriori that uncertainties in the gain 
calibration do not significantly affect our results. 

We argue that this may be due to relevant masses of ICM 
being pushed away perpendicularly to the direction of the 
merger in the trail of the bullet. Clearly, at this stage, we 


cannot discriminate from a full rotation due to an off-center 
merging from a disordered velocity field. We stress that if 
confirmed, this will be the detection of true, merger-induced 
bulk motion in a virialized cluster, at variance with a veloc¬ 
ity difference due to a pre-merger stage with two dark matter 
halos, still carrying their ICM, falling toward each other. 

In this work, we push the limits of what is possible with 
Chandra. In addition, the Bullet cluster constitutes a diffi¬ 
cult test for our method due to the very high temperatures 
involved, which make the measurement of the iron emission¬ 
line complex much harder than in colder ICM. Despite this, 
we find interesting hints as to the possible presence of bulk 
motions in the ICM due to the ongoing merger. Our future 
plan is to extend this study to a small sample of merging 
or post-merger clusters to identify the most promising cases 
where bulk motions can be successfully identified with Chan¬ 
dra. Currently, the capability of this analysis is limited due 
to the low spectral resolution of CCD data. In the future, the 
systematic measurement of the ICM velocity structure may 
have a strong impact on cluster physics, helping to reconstruct 
the assembly process of clusters across cosmic epochs and to 
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Fig. 15.— Best-fit energies for the fluorescent lines of Au La (upper panel) and Ni Ka (lower panel) lines in each region with the corresponding Icr error bars. 
Note that we are not able to identify the Ni Ka line in several regions, while the Au La line is always detected. 


derive more accurate total masses for cosmological tests. It 
would be important to plan future missions able to deliver spa¬ 
tially resolved observations with X-ray micro calorimeters to 
achieve a major breakthrough in the study of the dynamics of 
the ICM. 
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